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Abstract. The free-boundary compressible 1-D Euler equations with moving physical 
vacuum boundary are a system of hyperbolic conservation laws which are both charac- 
teristic and degenerate. The physical vacuum singularity (or rate-of-degeneracy) requires 
the sound speed c = 7p'^~^ to scale as the square-root of the distance to the vacuum 
boundary, and has attracted a great deal of attention in recent years [l]. We estab- 
lish the existence of unique solutions to this system on a short time-interval, which are 
smooth (in Sobolev spaces) all the way to the moving boundary. The proof is founded 
on a new higher-order Hardy- type inequality in conjunction with an approximation of 
the Euler equations consisting of a particular degenerate parabolic regularization. Our 
regular solutions can be viewed as degenerate viscosity solutions. 



Contents 

1. Introduction 1 

2. Notation and Weighted Spaces £ 

3. A higher-order Hardy-type inequality £ 

4. A degenerate parabolic approximation of the compressible Euler equations in vacuum S 

5. Solving the parabolic K-problem (|4.1|l by a fixed-point method S 

6. Asymptotic estimates for which are independent of 14 

7. Proof of Theorem O 24 

8. The case 7 7^ 2 26 
References 26 



1. Introduction 

1.1. Compressible Euler equations and the physical vacuum boundary. This paper 
is concerned with the evolving vacuum state in a compressible gas flow. For < t < T, 
the evolution of a one-dimensional compressible gas moving inside of a dynamic vacuum 
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in lit), 


(1.1a) 


in I{t), 


(1.1b) 


on r{t), 


(1.1c) 




(l.ld) 


on /(O), 


(l.le) 




(l.lf) 



boundary is modeled by the one- phase compressible Euler equations: 
p [ut + uu^] + p{p)x = 
Pt + {pu)x = 
p = 

v(r(t)) = u 

{p,u) = {po,uq) 
/(O) = / = {0 < x < 1} , 

The open, bounded interval I{t) C M denotes the changing domain occupied by the gas, 
r(t) := dl{t) denotes the moving vacuum boundary, V(r{t)) denotes the velocity of T{t). 
The scalar field u denotes the Eulerian velocity field, p denotes the pressure function, and 
p denotes the density of the gas. The equation of state p{p) is given by 

p{x,t) = p{x,t)"' for 7>1, (1.2) 

where Cj is the adiabatic constant which we set to one, and 

p > in I{t) and p = on r{t) . 

Equation (|l.lap is the conservation of momentum; (jl.lbp is the conservation of mass; the 
boundary condition (|l.lc|) states that the pressure (and hence density) must vanish along 
the vacuum boundary; (jl.ldp states that the vacuum boundary is moving with the fluid 
velocity, and (|l.lep - (ll.lf|) are the initial conditions for the density, velocity, and domain. 
Using the equation of state (|1.2p . p. lap is written as 

p[ut + uu^] + {p'>)^ in I(t) . dni) 
With the sound speed given by c^(a;, t) = ^p'^'^^{x, t), and with co — c{x, 0), the condition 



< 



dc. 



dx 



< oo on r (1.3) 



defines a physical vacuum boundary (or singularity) (see [TT] , [T^] , [13] , [H] ) ■ Since po > 
in /, (|1.3p implies that for some positive constant C and x G I near the vacuum boundary 

r, 

p^^'\x)>Cdist{x,T). (1.4) 
Equivalently, the physical vacuum condition (|1.4p implies that for some a > 0, 



-ix) 



> 1 for any x satisfying d{x, dl) < a , (1-5) 



dx 

and for a constant Cq, depending on a, 

Pq~^{x) > Ca > for any x such that d{x, dl) > a . (1.6) 

Because of condition (jl.4p . the compressible Euler system (II. ip is a degenerate and charac- 
teristic hyperbolic system to which standard methods of symmetric hyperbolic conservation 
laws cannot be applied in standard Sobolev spaces. In [2], we established a priori estimates 
for the multi-D compressible Euler equations with physical vacuum boundary. 

The main result of this paper is the construction of unique solutions in the 1-D case, 
which are smooth all the way to the moving vacuum boundary on a (short) time-interval 
[0,T], where T depends on the initial data. We combine the methodology of our a priori 
estimates [5], with a particular degenerate parabolic regularization of Euler, which follows 
our methodology in [31 H], as well as a new higher-order Hardy- type inequality which permits 
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the construction of solutions to this degenerate parabohc regularization. As we describe 
below in Section 11.71 our solutions can be thought of as degenerate viscosity solutions. The 
multi-D existence theory is treated in [5]. 

1.2. Fixing the domain and the Lagrangian variables on the reference interval 

/. Wc transform the system (jl.ip into Lagrangian variables. We let ri{x,t) denote the 
"position" of the gas particle x at time t. Thus, 

dti] ^ u o rj for t > and r/(a;, 0) = a; 

where o denotes composition so that [u o 7f\ {x, t) := u{ri{x, t),t) . We set 

V ^ uo rj (Lagrangian velocity), 

f — po (Lagrangian density) . 

The Lagrangian version of equations (|l.la|) - (jl.lbp can be written on the fixed reference 
domain / as 

M + (r)x = in/x(0,r], (1.7a) 

ft + fv,/Vx = in/x(0,r], (1.7b) 

/ = in/x(0,r], (1.7c) 

if,^,?]) ^ ipo,uo,e) in/x{< = 0}, (1.7d) 

where e(x) = a: denotes the identity map on /. 
It follows from solving the equation (jl.7bp that 

f ^ P°V = Po/Vx, (1-8) 

so that the initial density function po can be viewed as a parameter in the Euler equations. 
Let r := dl denote the initial vacuum boundary; we write the compressible Euler equations 
(EZl) as 

povt + {p2/v^x)x =0 in / X (0, T] , (1.9a) 

(?7,w) = (e,uo) in/x{t = 0}, (1.9b) 
p2~^^0 onL, (1.9c) 

with p2^^{x) > Cdist(x, F) for a; G / near F. 

1.3. Setting 7 = 2. We will begin our analysis of (|1.9p by considering the case that 7 = 2, 
in which case, we seek solutions ij to the following system: 

PoVt + {pl/riDx =0 in / X (0, T] , (1.10a) 

{r],v) ^ {e,uo) on/x{t = 0}, (1.10b) 
po=0 onF, (1.10c) 

with po{x) > Cdist(a;,F) for x € I near F. 
The equation (|1.10ap is equivalent to 

vt + 2rj-\poT]-'),^0, (1.11) 

and (jl.lip can be written as 

Vt+Po(.Vx% + nPo)xVx^ = 0. (1.12) 

Because of the degeneracy caused by po = on F, all three equivalent forms of the com- 
pressible Euler equations are crucially used in our analysis. The equation p.lOap is used 
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for energy estimates, while (jl.lip and (I1.12p are used for additional elliptic-type estimates 
that rely on our higher-order Hardy-type inequality. 

1.4. The reference domain /. As we have already noted above, the initial domain / C M 
at time t = is given by 

/=(0,1), 

and the initial boundary points are given hy T = dl = {0, 1}. 

1.5. The higher-order energy function for the case 7 = 2. We will consider the 
following higher-order energy function: 

4 2 
i?(t,z;) =^||aMt,-)ll^.-t(,)+^|lpo9t%(i,-)ll?f3-.(,) 

s=0 s=0 

+ W^dtdMt, -nhn) + WVf^d^dMt, ■ (1-13) 

We define the polynomial function Mq by 

Mo^ P{E{0,v)). (1.14) 

1.6. The Main Result. 

Theorem 1.1 (Existence and uniqueness for the case 7 = 2). Given initial data (uqiPo) 
such that Mq < 00 and the physical vacuum condition holds for po, there exists a 

solution to il.lU\) (and hence H-l]) ) on [0, T] for T > taken sufficiently small, such that 

sup E{t) < 2Mo . 
te[o,T] 

Moreover if the initial data satisfies 

3 3 

J2 wdMo, -m-s + E wpod^o, out. < 00 , (1.15) 

s=0 s=0 

then the solution is unique. 

Remark 1.2. The case of arbitrary 7 > 1 is treated in Theorem \8.1\ below. 

Remark 1.3. Given the regularity provided by the energy function hl.l3\) . we see that the 
Lagrangian flow map rj € C([0, T], H^{I)). In our estimates for the multi-D problem [2], we 
showed that r/ gains regularity with respect to the velocity field v, and this fact is essential to 
control the geometry of the evolving free-surface. This improved regularity for rj also holds 
in the TD setting, but it is not necessary for our estimates, as no geometry is involved. 

Remark 1.4. Given uq and po, and using the fact that 77(2;, 0) = x, the quantity vt\t=o is 
computed using 1 1.10a]) : 



po 



= -2— 

t=o dx 



Similarly, for all k £ 



aM,=o = - ( -{plhl' 

\Pa 



so that each dfV\t=o is a function of space- derivates of uq and pQ. 
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Remark 1.5. Notice that our functional framework provides solutions which have op- 
timal Sobolev regularity all the way to the boundary. Hence, in the physical case that 
c ^ ^ dist{dl), no singular behavior occurs near the vacuum boundary, even though both 
families of characteristics cross, and in particular, meet tangentially to T{t) at a point. 

Remark 1.6. Because of the degeneracy of the density function po at the initial boundary 
dl , no compatibility conditions are required of the initial data. 

1.7. History of prior results for the compressible Euler equations with vacuum 
boundary. Some of the early developments in the theory of vacuum states for compressible 
gas dynamics can be found in [14j and [9]. We are aware of only a handful of previous 
theorems pertaining to the existence of solutions to the compressible and undamped Euler 
equations with moving vacuum boundary. Makino [15j considered compactly supported 
initial data, and treated the compressible Euler equations for a gas as being set on M'^ x (0, T] . 
With his methodology, it is not possible to track the location of the vacuum boundary (nor is 
it necessary); nevertheless, an existence theory was developed in this context, by a variable 
change that permitted the standard theory of symmetric hyperbolic systems to be employed. 
Unfortunately, the constraints on the data are too severe to allow for the evolution of the 
physical vacuum boundary. 

In [10], Lindblad proved existence and uniqueness for the 3-D compressible Euler equa- 
tions modeling a liquid rather than a gas. For a compressible liquid, the density p > A > is 
assumed to be a strictly positive constant on the moving vacuum boundary T{t) and is thus 
uniformly bounded below by a positive constant. As such, the compressible liquid provides 
a uniformly hyperbolic, but characteristic, system. Lindblad used Lagrangian variables 
combined with Nash-Moser iteration to construct solutions. More recently, Trakhinin [T7] 
provided an alternative proof for the existence of a compressible liquid, employing a solution 
strategy based on symmetric hyperbolic systems combined with Nash-Moser iteration. 

In the presence of damping, and with mild singularity, some existence results of smooth 
solutions are available, based on the adaptation of the theory of symmetric hyperbolic 
systems. In [T^], a local existence theory was developed for the case that c" (with < a < 1) 
is smooth across F, using methods that are not applicable to the local existence theory for 
the physical vacuum boundary. An existence theory for the small perturbation of a planar 
wave was developed in [18] . See also [13| and , for other features of the vacuum state 
problem. 

The only existence theory for the physical vacuum boundary condition that we know of 
can be found in the recent paper by Jang and Masmoudi 7J for the 1-D compressible gas, 
wherein weighted Sobolev norms are used for the energy estimates. From these weighted 
norms, the regularity of the solutions cannot be directly determined. Letting d denote the 
distance function to the boundary dl, and letting || • ||o denote the L^(/)-norm, an example 
of the type of bound that is proved for the velocity field in JT] is the following: 

\\dv\\l + ||rfWa;||o + \\dVxx + 2vx\\l + \\dVxxx + '2Vxx - 2d"^ VxWl 

<oo. (1.16) 

The problem with inferring the regularity of v from this bound can already be seen at the 
level of an H^{I) estimate. In particular, the bound on the norm + 2Da;|lg only implies 

a bound on ||du2;2;||o and ||wa;|jg if the integration by parts on the cross-term, 



4 J dvxxVx = 



-2j^dx\Vx\^, 
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can be justified, which in turn requires having better regularity for Vx than the a priori 
bounds provide. Any methodology which seeks regularity in (unweighted) Sobolev spaces for 
solutions must contend with this type of issue. We overcome this difficulty by constructing 
(sufficiently) smooth solutions to a degenerate parabolic regularization, and thus the sort 
of integration-by-parts difhculty just described is overcome. One can view our solutions as 
degenerate viscosity solutions. The key to their construction is our higher-order Hardy-type 
inequality that we provide below. 

2. Notation and Weighted Spaces 

2.1. Differentiation and norms in the open interval /. Throughout the paper the 
symbol D will be used to denote 

For integers fc > 0, we define the Sobolev space H'^{r} to be the completion of C°°(/) in 
the norm 

1/2 



I \D''u{x)\^dx 



For real numbers s > 0, the Sobolev spaces H'^{I) and the norms || • ||s are defined by 
interpolation. We use Hq{I) to denote the subspace of H^{I) consisting of those functions 
u(x) that vanish at a; = and x ~ I. 

2.2. The embedding of a weighted Sobolev space. Using d to denote the distance 
function to the boundary F, and letting p = 1 or 2, the weighted Sobolev space H^p (/), with 
norm given by Jj d{x)P {\F{x)\'^ + \DF{x)\'^) dx for any F e H\p{I), satisfies the following 
embedding: 

so that there is a constant C > depending only on / and p such that 

ll^ll?-p/2 < C jd{xY{\F{x)\' + \DF[x)\') dx. (2.1) 
See, for example, Section 8.8 in Kufner [5]. 

3. A HIGHER-ORDER HARDY-TYPE INEQUALITY 

We will make fundamental use of the following generalization of the well-known Hardy 
inequality to higher-order derivatives: 

Lemma 3.1 (Higher-order Hardy-type inequality). Let s > 1 be a given integer, and suppose 
that 

ue H'{i)nH^{i). 

Then if d denotes the distance fuction to dl , we have that - e H^^^ll) with 

d 



<C\\u\\,. (3.1) 

-1 

Proof. We use an induction argument. The case s = 1 is of course the classical Hardy 
inequality. Let us now assume that the inequality (j3.ip holds for a given s > 1, and suppose 
that 

u^H'+^{I)r^Hl{I). 

Using D to denote a simple computation shows that for to g N, 

^'"(^) = ^- (3.2) 
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with 

m 



k=0 



for a constant depending on k and m. From the regularity of u, we see that / G Hq{I). 



Next, with D = we obtain the identity 



fc=0 fe=0 



s-l 



+ C^+^D'~''u {k + l)\{~l fd'^''-^ 

k=0 

= D'+^u s\{~lYd' . (3.3) 
Since / G Hq{T), we deduce from p.3p that for any x e (0, 

/(x) = (-l)^s! r D^+'uiy)y^dy, 



which by substitution in (13. 2|) yields the identity 



which by a simple majoration provides the bound 

rf ' d(a;) 

where tpi is the piecewise affine function equal to 1 on [0, and to on [|, 1]. Next, for 
any x € [55 l)i we obtain similarly that 

d d{x) 

where ■02 is the piecewise affine function equal to on [0, i] and to 1 on [|, 1]. so that for 
any x & I: 

lr.sr^,r M ^ I M^) lo \D'+My)\ dy + M^)I^\D'+My)\ dy 
\u{^}{x)\ss. ^^^^ . 

Now, with g = \D''+'^u{y)\ dy + ^'2(2:) \D''+'^u{y)\ dy, we notice that g e H^il), 

with 

\\g\\i<C\\D^+'u\\o. 

Therefore, by the classical Hardy inequality, we infer from p.4p that 

\\D%^)\l<C\\gh<C\\D^+'u\\o. (3.5) 

Since we assumed in our induction process that our generalized Hardy inequality is true at 
order s, we then have that 

\OU<C\\u\\s, 
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which, together with (13. 5p . impHes that 



and thus estabhshes the property at order s + 1, and concludes the proof. □ 

In order to obtain estimates independent of a regularization parameter k, defined in 
Section [4l we will also need the following Lemma, whose proof can be found in Lemma 
1, Section 6 of [3]: 

Lemma 3.2. Let k > and g e L°°(0, T; ''(/))) be given, and let f £ ^^(0, T; be 
such that 

f + Kft^g in{0,T)xI. 

Then, 

ll/l|L~(o,T;_f/=(7)) <C max{||/(0)||s, ||.9||l=o(o,T;//=(7))}- 

4. A DEGENERATE PARABOLIC APPROXIMATION OF THE COMPRESSIBLE EULER 

EQUATIONS IN VACUUM 

For K > 0, we consider the following nonlinear degenerate parabolic approximation of the 
compressible Euler system (|1.10l) : 

Povt + {plhl)x = h[pIvx\x in / X (0, T] , (4.1a) 
(?7,w) = (e,uo) on / x = 0} , (4.1b) 
Pq=0 onT, (4.1c) 

with po{x) > C dist(a;. F) for x G / near F. 

We will first obtain the existence of a solution to (|4.ip on a short time interval [0, T„] 
(with depending a priori on k). We will then perform energy estimates on this solution 
that will show that the time of existence, in fact, does not depend on k, and moreover 
that our a priori estimates for this sequence of solutions is also independently of k. The 
existence of a solution to the compressible Euler equations (|1.10p is obtained as the weak 
limit as K — > of the sequence of solutions to (|4.ip . 

5. Solving the parabolic k-problem (14. 1|) by a fixed-point method 
For notational convenience, we will write 



/ dri 



and similarly for other functions. 



5.1. Assumption on initial data. Given uo and po, and using the fact that rj{x,Q) — x, 
the quantity vt\t=o for the degenerate parabolic K-problem is computed using (|4.1a|) : 



\Pa Po 



-[plu',r-2p', 
Po 



Similarly, for all fc e N, 



t=0 



These formulae make it clear that each dfV\t=o is a function of space-derivates of uq and po- 



MOVING-BOUNDARY 1-D COMPRESSIBLE EULER EQUATIONS IN VACUUM 



9 



For the purposes of constructing solutions to the degenerate parabolic K-problem ()4.ip . 
it is convenient to smooth the initial data. By standard convolution methods, we assume 
that uo g C°°{I) and that po G C°°{I) and satisfies (HHD and (fTB)) . 

5.2. Functional framework for the fixed-point scheme and some notational con- 
ventions. For T > 0, we shall denote by Xt the following Hilbcrt space: 

XT = {ve L\Q,T-H\I))\ dfv e L^{0,T;H\I)), podfv G L\Q,T-H\I)) (5.1) 

dlveL\Q,T-H\l)), pod!veL\0,T;H'{I))}, (5.2) 

endowed with its natural Hilbcrt norm: 

M\l2(o,T;HHi)) + \\9tv\\l2(o.T:HHi)) + \\Podfv\\] 



U^IIXt — ll'-'llL2(0,T;H2(/)) -t- \\Uf t'||L2(0,T;ffi(/)) llW^t ^llL2(0,T;if2(/)) 
+ \\dfv\\L2(QT-H^(I)) + \\Podfv\\l2(oT-H='iI)) ■ 

For M > given sufficiently large, we define the following closed, bounded, convex subset 
of A^r: 

Ct(M) = {v€Xt : \\v\\x^ < M}, (5.3) 

which is indeed non empty if M is large enough. Henceforth, we assume that T > is given 
such that independently of the choice of t; G Ct{M), 

ri{x,t)=x+ / v{x,s)ds 
Jq 

is injective for t G [0,T], and that ^ < ■q'{x,t) < | for t G [0,r] and a; G /. This can be 
achieved by taking T > sufficiently small: with e{x) = x, notice that 



||7/(-,i)-e||i = II / v'{-,s)ds\\^ < VtM. 
Jo 



The space Xt will be appropriate for our fixed-point methodology to prove existence of a 
solution to our ^-regularized parabolic problem (|4.ip . 

Finally, we define the polynomial function A/q of norms of the initial data as follows: 

Mo = PiW^dfv'mio, ||a^(o)||i, ||a^(o)||2, ||9^(o)||2, WdtvmU, hoh) . (5.4) 

5.3. A theorem for the existence and uniqueness of solutions to the parabolic 

K-problem. We will make use of the Tychonoff fixed-point Theorem in our fixed-point 
procedure (see, for example, ^6^). Recall that this states that for a reflexive separable Banach 
space Xt, and Ct{M) C Xt a closed, convex, bounded subset, if F : Ct{M) Ct{M) is 
weakly sequentially continuous into Xt, then F has at least one fixed point. 

Theorem 5.1 (Solutions to the parabolic K-problem). Given our smooth data, forT taken 
sufficiently small, there exists a unique solution v G Xt to the degenerate parabolic K-problem 

5.4. Linearizing the degenerate parabolic K-problem. Given v G Ct{M), and defining 
f]{x, t) ^ X + v{x, r jdr, we consider the linear equation for v: 



Povt 



77'2 



K[piv']'^Q. (5.5) 



We will prove the following: 

(1) V \s a unique solution to l|5.5p : 

(2) V G Ct{M) for T taken sufficiently small; 
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(3) the map v t—^ v : CxiM) — + Ct{M), and is sequentially weakly continuous in Xt- 
The solution to our parabolic K-problem (|4.ip will then be obtained as a fixed-point of the 
map w I— > w in Xt via the Tychonoff fixed-point Theorem. 

In order to use our higher-order Hardy-type inequality, Lemma |3.1[ it will be convenient 
to introduce the new variable 

X = pov' 

which then belongs to H^{I). By a simple computation, we see that (jS.Sp is equivalent to 



"2 /poV' 


1 


■ 1 




r 




— K 








f]'\fj' ) 




.Po 







= 0, 



and hence that 



PO 



-(po X)' 

PO 



in [0, T\ X /, 

on [0, T] X a/, 
on /. 



(5.6a) 

(5.6b) 
(5.6c) 



X = 

We shall therefore solve the degenerate linear parabolic problem ()5.6|) with Dirichlet bound- 
ary conditions, which (as we will prove) will surprising admit a solution with arbitrarily high 
space regularity (depending on the regularity of the right-hand side and the initial data, of 
course), and not just an i7Q(/)-type weak solution. After we obtain the solution X, we will 
then easily find our solution v to ()5.5|) . 

In order to construct our fixed-point, we will need to obtain estimates for X (and hence v) 
with high space regularity; in particular, we will need to study the fifth time-differentiated 
problem. For this purpose, it is convenient to define the new variable 

Y = dlX = p^d\v' , 

and consider the following equation for Y: 



(5.7) 



Po 



-(po Y)' 

Po 



Y\ 



Y = 

=0 — ^init 



in [0, T] X /, 

on [0, T] X dl, 
in / , 



(5.8a) 

(5.8b) 
(5.8c) 



where Y„ 



Podfv'lt^Q. 



5.5. Existence of a vireak solution to the linear problem (|5.8p by a Galerkin 
scheme. Let (e„)„gN denote a Hilbert basis of Hq{I), with each e„ being of class H^{I). 
Such a choice of basis is indeed possible as we can take for instance the eigenfunctions of 
the Laplace operator on / with vanishing Dirichlet boundary conditions. We then define the 
Galerkin approximation at order n > 1 of (|5.6p as being under the form F„ — X]r=o 
such that: 



Vfce{0,...,fc}, ( ^ + K(por„)' , ^ 
Po Po 



(i^nit, ei)L2(/). 



in[0,T], (5.9a) 
(5.9b) 



Since each is in H^{I) n Hq{I), we have by our high-order Hardy- type inequality (|3.1 
that 

Po 
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therefore, each integral written in (15. 9p is well-defined. Furthermore, as the Ci are linearly 
independent, so are the and therefore the determinant of the matrix 



(( 



3_\ 



/PO VPO ^'V(»,i)GN„={l,...,n} 

is nonzero. This implies that our finite-dimensional Galerkin approximation (|5.9[) is a well- 
defined first-order differential system of order n + 1, which therefore has a solution on a time 
interval [0, T„] , where T„ a priori depends on the dimension n of the Galerkin approximation. 
In order to prove that r„ = T, with T independent of n, we notice that since Yn is a linear 
combination of the (i G N„), we have that 

- n[^{po YJ]', Y^ = (dn^[^]'] , y,:) . (5.10) 
\Po Pa J LHi) V ^V' V' J L^i) 

Since Yn G Hq{I) and [^(po ^n)']' ^ -^^(-^); integration by parts yields 

-(po i;.)']' Yn^ [ [-{po r„)']K = / + I p'o-Y:,. (5.11) 
/ po J/ po J/ J I po 

Next, using our higher-order Hardy-type inequality, we see that ^ e H^{I), and thus 



/ P'o-YL f p'o-YX. + I '\yI - /po- 
h Po Ji Po J I Po J I Po 

which implies that 

j,p'pJ-~ 2], pf- 2i/° 

Substitution of this identity into (|5.10p and (|5.1ip yields 

^ ri r Y'^ r Y'^ r \ 1 0'"^ 



II ^ n 



which shows that (since our given v e Ct{M)): 



d r Y"^ r Y^ r r y"^ 



J I Po J I Po J I J I Po 

for a constant C„ depending on k. Consequently, T„ = T with T independent of n, and 

sup / ^ + « / / < c-.T -f CAAo , (5.12) 

[0,T] J/ Po Jo J I 

where No is defined in (|5.4p . Thus, there exists a subsequence of (y„) which converges 
weakly to some Y E L^(0, T; (/)), which satisfies 

sup/ — + k/ [ Y'^ < C^T + CMo . (5.13) 

[0,T] J/ Po Jo J I 

With (|57ll . we see that 

\\pod!v'\\L2io,T;mii)) < C^T + CN-Q. (5.14) 
Furthermore, it can also be shown using standard arguments that F is a solution of (|5.8p 
(where (|5.8ap is satisfied almost everywhere in [0, T] x / and holds in a variational sense for 
all test functions in L^(0, T; iJp (/))), and that 

^eL^Q,T;H-\l)). 
Po 
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Now, with the functions 



^ _ d'v' 



for i = 0,1,2,3,4, 



t=o 



we define 



and 



X{t,x)=J2^^*' 



Jo Jo Jo 



Z{ti, x)dtidt2dt^dt^ 



(5.15) 



(5.16) 



We then see that X e C"([0, T];H^{I)) is a sohition of jSj]), with d^X = Y. In order to 
obtain a fixed-point for the map w i-^ u, we need to establish better space regularity for Z, 
and hence X and v. 

5.6. Improved space regularity for Z. We introduce the variable v defined by 

X{t,-) 



Jo Po(-) 

so that V vanishes at x = 0, and wiU us to employ the Poincare inequality with this variable. 
It is easy to see that 

X = pqv', 

and 

Z = podfv'. 

Thanks to the standard Hardy inequality, we thus have that 

liayiio<c||z||i, 

and hence by Poincare's inequality, 



(5.17) 



With 



\\dtv{t,-)\\,<C\\Z{t,-)\\ 



po rj' ij 



(5.18) 



our starting point is the equation 

- - 4- iPo Zyy = -df [4 [?] '] '+Fo in [0, T] x /, 
Po Po n ^ri" ' 

which follows from our definition of Z given in (j5.15p and time- integration of (|5.8ap . From 
this equation, we infer that 

4[-ip^ ^)rii„<||5,n4[?]Tllo + ll-llo + ii^oiio. 

Po L^i L ^1 J J \n) " pq"o 

By the standard Hardy inequality and the fact that v e Ct{M), we obtain the estimate 

k\\[-{po zyfW < CMVT + C\\Y\\i+Afo, 
Po ^ 

where Cm is a constant that depends on M. In particular, using (|5.17|) . we see that 

-ipo Zy^podtv" + 2p',dfv' 
Po 

so that 



K\\podtv"' + 3p'odtv" + 2p'idtv'\l < CmVt + CWYWi+Mo , 
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which impUes that 

«||(po5tV'||o < CmVt + C\\Y\\^ + AAo + n\\p'^'dtv\\o + KW^dtv'W^ 

< CmVt + C\\Y\\i +Mo + <\\p'o\\l^ + 3||p^'|Uoe)||Z||i , 

where we have used (15.18^ in the second inequaUty above. Having estabhshed in (|5.18p that 
{podfv) £ Hq{I), eUiptic regularity shows that 



4Podfvh < CmVt + C\\Y\\i + AAo + aj(||p'o"||loo + 3||p;,'||l==)||^||i ■ (5.20) 
Now, thanks to our high-order Hardy-type inequahty, we infer from fr'om (|5.20p that 

^dfvh < CmVt + C\\Y\\i +Afo + «(|1po"IU- + 3||Polk-)ll^lli • (5-21) 
Next we see that (|5.20p implies that 

K\\podtv' + p'dtvh<CMVT + C\\Y\\i+Afo + K{\\p'^'\\L^+3\\p'^\\L^)\\Z\^^ 
which thanks to (|5.2ip and (|5.17p implies that 

4Z\\2 < CmVt + C\\Y\\i + Mo + h{\\p'^'\\l^ + 3\\p'^\\l^)\\Z\\i . (5.22) 

5.7. Definition of v and the existence of a fixed-point. We are now in a position to 
define v in the following fashion: let us first define on [0, T] 



fit) = uo{0) - 



r-[4]'(-'0)+«^ f-[poxn;0), 

la Po V Jq po 

which is well-defined thanks to (|5.22p and (|5.13p . We next define 

v{t,x) ^ f{t)+v{t,x). (5.23) 
We then notice that from (|5.6p . we immediately have that 

from which we infer that in [0, T] x I 



Po^v''^' Po' 



for some function g depending only on t. By taking the trace of this equation on the left 
end-point a; = 0, we see that 

vtim + l[4]'(i,0) - K-[plv'nt,0) = g{t), 



Po "V'^ ■ Po' 



which together with the identity 

vt{t,0) = Mt) = -l[^]'{t,0) + K-[plv'nt,0) 
Po V Po 

shows that g{t) = 0. Therefore, w is a solution of (|5.5p . and also satisfies by construction 
viO,-)=UQ{-). 

We can now establish the existence of a fixed-point for the mapping u — s- t; in Ct{M), 
with T taken sufficiently small and depending a priori on k. We first notice that, thanks to 
the estimates (|5.22p and (|5.13p . we have the inequality 
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which together with ()5.2H) and (I5.13P provides the estimate 

I|5M|l2(o,t;H=(7)) < AAo + C.VtCm . (5.24) 

Then, (|O0)) imphes that 

||aMU.(o,T;H2(/)) < AAo + C.VTCm , (5.25) 
and combining this with (|5.14p shows that 

\\Pod!v'\\mo.T,mii)) < AAo + C^VtCm , (5.26) 

In turn, (|5.25p shows that for 

r<7f7^, (5.27) 

V e Ct(M). Moreover, it is clear that there is only one solution v e L^{0,T; H^{I)) of 
(|5.5p with z;(0) = uq (where this initial condition is well-defined as ||wt||i2(o,T;iri(/)) 5: 
Ao + CkVTCm), since if we denote by w another solution with the same regularity, the 
difference Sv = v — w satisfies dv{0, = with poSvt — k[pqSv']' = which implies 

1 d f 2 f 2 r 2 

2dtJ/'^'' +-J/oSv =0, 

which with 5v{Q, •) = implies Sv — 0. So the mapping u ^ u is well defined, and thanks to 
(j5.24p is a mapping from Ct{M) into itself for T — T^ satisfving |5.27|) . As it is furthermore 
clear that it is weakly continuous in the L'^{0,Ti^; H^{I)) norm, the Tychonoff fixed-point 
theorem [6] provides us with the existence of a fixed-point to this mapping. Such a fixed- 
point, which we denote by w^, is a solution of the nonlinear degenerate parabolic /c-problem 
(|4.ip . with initial condition Wk(0, •) = wo(')- K should be clear that the fixed-point also 
satisfies (|5.14p and (|5.25p so that 

l|Po9fuK||L2(o,T://3(/)) + l|3t^«K||L2(o,T://2(/)) 

+ \\pod^v'jL2(^o,T-HHI)) + \\d!v,\\L2^o^T;HHI)) <M, (5.28) 

In the next section, we establish ^-independent estimates for in L^{0, T^,; H'^{I)) (which 
are indeed possible because our parabolic approximate K-problem respects the stucture of 
the original compressible Euler equations (ll.lOp ). from which we infer a short time- interval 
of existence [0,T], with T independent of k. These ^-independent estimates will allows us 
to pass to the weak limit of the sequence as k ^ to obtain the solution to (jl.lOp . 

6. Asymptotic estimates for u„ which are independent of k. 

6.1. The higher-order energy function appropriate for the asymptotic estimates 

as K — > 0. Our objective in this section is to show that the higher-order energy function E 
defined in (|1.13p satisfies the inequality 

sup E{t) < Mq + CT P{ sup E{t)), (6.1) 
te[o,T] te[o.T] 

where P denotes a polynomial function, and for T > taken sufficiently small, with Mq 
defined in p.l4p . The norms in E are for solutions v^, to our degenerate parabolic K-problem 
()4.ip . According to Theorem l5.H G Xt^ with the additional bound ||9('uK||L2(o,T;^;if2(/)) < 
oo provided by (|5.24p . As such, the energy function E is continuous with respect to t, and 
the inequality (|6.ip would thus establish a time interval of existence and bound which are 
both independent of n. For the sake of notational convenience, we shall denote u„ by v. 
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6.2. A K- independent energy estimate on the fifth time-differentiated problem. 

Our starting point shall be the fifth timc-diffcrcntiated problem of (14. 1|) for which we have, 
by naturally using dfv £ i^(0, T^; H^{I)) as a test function, the following identity: 



^dt Jj 



i5„-v|2 



Po\d!v 



(6.2) 



l2 



In order to form the exact time-derivative in term Ti , we rely on the fact that solutions we 
constructed to (|4.ip satisfy ^^v £ L^(0, T„; which follows from the relation 





'2 /poy 


K 

H 






77' V77'/ 


PQ 





and the estimate (j5.28p . Upon integration in time, both the terms Xi and X3 provide sign- 
definite energy contributions, so we focus our attention on the nonlinear estimates required 
of the term l2- 
We see that 

-I. = 2 / dfv'[^]d!v' + Y.<^^ f Ot^ldtUdfv' 
Ji I Ji I I 



d 

dt J J 



Hence integrating (|6.2p from to t £ [0, T^], we find that 



„2/q5 -/\2 



Jl 

^i~i\2 Pa 



T] 



JI 



a=l Jl V V 



(6.3) 



We next show that all of the error terms, comprising the right-hand side of (|6.3p can be 
bounded by Ct P (snp E). For the first spacetime integral appearing on the right-hand 
side of (16. 31). it is clear that 



3 f f{dfvrv'^<CtP{snpE). 
Jo Ji V [OA] 



(6.4) 



We now study the last integrals on the right-hand side of (|6.3p . Using integration- by-parts 
in time, we have that 



lo Ii 



Jo Ji V V J I V V 



(6.5) 



We first consider the spacetime integral on the right-hand side of (|6.5p . As the sum is taken 
for a = 1 to 4, we then see that it will be written under the form of the sum of spacetime 
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integrals of the following type: 



h - J^PodfvR{fi)podtv', 

h^J^ J^Pod!v'v'dti'R{fi)podtv', 

h = f J^pod?i'v'd^v'R{fi)podti\ 
t p 

!2~l I r>/~\/Ci ~!\2 a4~? 



/4 = j^podtv'v'Rmdtv'Yp^dtv' , 
Where R{fj) denotes a rational function of fj. We first immediately see that 



|/i I < Ct P(supi;). (6.6) 

[o,t] 



Next, we have that 



Similarly, 



\h\ <c j'^ jy^dlv'\\LA\v^\L^\\dtv'\\LA\R{mL^\\p^dtv 

<C fju,d^i'\\^.\\v'\\L^\\d,i'\\^.jR{mL^^^^^ 

<CtP{supE). (6.7) 

[0,t] 



|/3| <C /Voa,V|U^||z;'|U^||a,V||o||i?(77)|k~||poa,V||o 
Jo 

<CtP{supE), (6.8) 

[o.t] 



and 



\h\ <C f ||po9*V|lL~|k'|lL~|ia,5'|li.||i?(77)|U~||po5,V|lo 
Jo 

<CtP{sn]yE), (6.9) 
[o,t] 

where we have used the fact that in 1-D, 11 • ||l°° < C\\ ■ ll^fi and 11 • 11^4 < CII • 11 i . Therefore, 
estimates (|6.3p - (|6.7p provide us with 

\ ( podfv'it) + [ idtv'r4(*) + I Plidti'f < Mo + Ct P{snpE), 
2 J I J I Jo J I [o,t] 

and thus, employing the fundamental theorem of calculus. 



podti'it) + ipodti'Yit) + n piidtv') 



Ji 



< Mo + Ct P{snpE) + 3 / [dfv'fWpl f 4 

[O.t] Ji Jo V 

< Mo + Ct P{snpE). (6.10) 

[o,t] 
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6.3. Elliptic and Hardy-type estimates for d^v{t). Having obtained the energy es- 
timate ()6.10p for the fifth time-differentiated problem, we can begin our bootstrapping 
argument. We now consider the third time-differentiated version of (|4.1ap . 



which can be written as 



and finally rewritten as the following identity: 

-2 [pld?v'] ' - n[pldfv']' = - podtv + c, [^^] ' + c, [^] ' 

-2[p^5,V]'(l-i,)-6,ga,V|;. (6.11) 

Here, ci and C2 are constants whose exact value is not important. 

Therefore, using Lemma 13.21 together with the fundamental theorem of calculus for the 
fourth term on the right-hand side of (|6.1ip . we obtain that for any t G [0,T„], 

II 2 r 2o2~n'|| ^ ||q4~|| , \\ ^1 r Pldtv' v' y . C2 rpgw'^ -, 

[o,t] Vo' ' io,t] [o,t] Vo' r " [0,t]"PO 

+ sup||-[p292^']'||o||3 / 4||i^+6sup||poatV^||Q . (6.12) 
[o,t]"Po Jo V [o,t]" V 

We next estimate each term on the right hand side of ()6.12p . For the first term, we will use 
our estimate (j6.10p from which we infer for each i G [0, Tf^]: 

[ plidfv^ + dtv'^]{t) < Mo + Ct P{supE). 

Note that the first term on the left-hand side of (|6.13p comes from the first term of (|6.10p . 
together with the fact that dfv{t, x) = V4{x) + / dfv{-, x). Therefore, the Sobolev weighted 

"'0 

embedding estimate ()2.ip provides us with the following estimate: 

/ dfi^it) < Mo + Ct P{supE). (6.13) 
Ji [o,t] 
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The remaining terms will be estimated by simply using the definition of the energy function 
E. For the second term on the right-hand side of ()6.12|) . we have that 

\\-[ fj^A \ Wo ^ IK^o^t) lloll^lLoc + ll^tL^J llo 

< CWipoVt) \\o\\v 111 + \\vt[^\ llo + ht [^\ llo + 4||«t llo 



<c\\{pov[y+ / {pov'uy\\o\\vT!H+ / vm + cH + / «;,iioii^'ii3 

Jo Jo ^ Jq * 

+ C\\v'M\Pov"U+C\\v'U\\ [ v"\\o\\pov[ + [ PoKtWi 



<\\{Pov[r+ / ipov',,y\\oW\\IH + 
Jo 

+ C\\v[ + [ ^llollS'llfll^^'i- 



+ CU\\o\\pov"\\t\\pov'^+ / Pov't'WL 

Jq 

where we have used the fact that || • \\l<x < C\\ ■ Iji. Thanks to the definition of E, the 
previous inequality provides us, for any t £ [0,T„], with 

supII — [^i^^l'|l < Csup£;3 {Mo + tP {sup E)). (6.14) 

[OA]" Po ' [0,i] [OA] 

For the third term on the right-hand side of (|6.12[) . we have similarly that 

ii^[f^]ii.<iita.')'iuHi|;L.+ii«'[^]'ii„ 



fjiS J llo II L fji5 J llo II L jy6 

<c\\{pov'„y+ I {pov[y\\h\ipov'y\\hK+ I y't\\l + cu + [ 







■C\\i'\\l\\Povl\L^+Cml\\Pov''\\L^ 



<c\\{pov'oy+ {pov'tyu\\{pov'y\\tH+ 5^111 + qi^, 







+ CH + I K\\l\\pov'^+ I pov','\\l\\pod"\\l+CH+ I v'MW I pov"\h, 
Jo ^ Jo Jo ^ Jo 

where we have used the fact that || • \\lp < C'p || • ||i, for all 1 < p < oo. Again, using the 

definition of E, the previous inequality provides us for any t € [0, T^] with 

1 o'^v'^ I 1 

sup|| — r^-^l II <Csupi;2 {Mn+tP{sviY>Ey). (6.15) 
[o,t] Po ?r " [OA] [o,t] 

For the fourth term of the right-hand side of (|6.12p . we see that 

\\^[pld^^i']'Ui [ y-\\^^{t)<C[ ||poa,V'||o + ||9y||o]tsup||5||2 
Po Jo V [o,t] 

<CtP{supE). (6.16) 

[OA] 
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Similarly, the fifth term of the right-hand side of (I6.12p yields the following estimate: 



<C\\pndti'U / v"h 
Jo 

< CtP{snpE) . (6.17) 

[o.t] 

Combining the estimates (|6.14p - (|6.17p . we obtain the inequality 

snp\\ — [pld^v']'L<CtP{supE) + CsupEi {Mq + tP{snp E)). (6.18) 

[0,t] PO [OA] [o,t] [o,t] 

At this stage, we remind the reader, that the solution v to our parabolic K-problem is in 
Xt^, so that for any t £ [0,Tf,], dfv £ H^{I). Notice that 

Po 

so (|6.18p is equivalent to 

supllpo^t^w" + 2p'od^v'\\ < CtP{supE) + CsupE^ (Mq + tP{supE)). (6.19) 

[0,t] [O.t] [O.t] [O.t] 

From this inequality, we would like to conclude that both ||9j'^t;'||o and ||poC^t^i'"||o are 
bounded by the right-hand side of (j6.19p : the regularity provided by solutions of the k- 
problem allow us to arrive at this conclusion. 

By expanding the left-hand side of (|6.19p . we see that 

sup||po9t'^" + 2p'od!v'\\l = \\podfv"\\l + Wodfi'Wl + 4 / Pod^i" p',d^,v' . (6.20) 
[o,t] Ji 

Given the regularity of d^v provide by our parabolic K-problem, we notice that the cross- 
term in (|6.20p is an exact derivative, 

:^2~," / £12--./ _ n I „ J la2„-./|2 



p,dtv" p'.dtv' = 2j^p,p',fjdtv'\\ 

so that by integrating-by-parts, we find that 

aJ^PoO^v" p'od^i' = -2\\p',d^i'\\l- J^pod^v' p'^d^v', 

and hence (|6.20p becomes 

sup||po5,^C" + 2p',dfv'\\l = Wpod^Wl + 2\\p'odfv'\\l - [ Pod^v' p'^d^v' . (6.21) 
[o.t] Ji 

Since the energy function E contains podfv{t) £ H^{I) and dfv{t) £ H^{I), the fundamental 
theorem of calculus shows that 

/ podfv' podfv' < CtP{supE) + CsupEi {Mq + tP{supE)) . 
Ji [o.t] [o,t] [o,t] 

Combining this inequality with (|6.2ip and (|6.18p . yields 

sup[|ipoc't^«"llo + l|Po<9t^*'llo] < CtP(sup£;) + Csupi;t {Mq + tP{snp E)), 

[0,t] [0,t] [0,t] [0,t] 
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and thus 

snp[\\pod^v"\\o + \\p'od^i'\\o + \\pod?i'\\o] < Mo + CtP{snpE) + CsupEi {Mo +tP (sup E)), 

[0,t] [0,t] [0,t] [0,t] 

and hence with the physical vacuum conditions on po given by (jl.Sp and (|1.6p . we have that 
sup[\\pQd'fi"\\o + Wdfi'Wo] < Mq + CtP{snpE) + CsiipEi (Mo + tP{supE)), 

[0,t] [0,t] [0,t] [0,t] 

which, together with (j6.13p . provide us with the estimate 

sup[\\pod^v"\\o + Wd^dWi] < Mo + CtP{supE) + CsupEi {Mo + tP{supE)). (6.22) 

[0,t] [0,t] [0,t] [0,t] 

6.4. Elliptic and Hardy-type estimates for v{t). Having obtained the estimates for 
dfv{t) in (|6.22p . we can next obtain our estimates for v{t). To do so, we consider the first 
time-difi^erentiated version of (|4.1ap . which yields the equation 

-2[^]'-4pld,v'r = -podfv, 

which we rewrite as the following identity: 

-^[pli']'- ^[pld^i'r = - dfi - A - 4) . (6.23) 
Po Pq Po V 

Using Lemma [3?^ we see that for any t G [0,Tk], 

snp\\-[pli']'\\^<Csnp\\dM, + Csnp\\-[pli']'U [ + C\\-[pli'r .^^ , 

m" Pa [o,t] [o,t] Po' Jo V Po 



(6.24) 

We next estimate each term on the right hand side of (|6.24p . The bound for the first term 
on the right-hand side of (|6.24p is provided by (|6.22p . 

The second term of the right-hand side of (|6.24p is estimated as follows: 

||-[Po^']'llill / i\\L^it)<C[\\pov"'\\o + \\vh]tsnp\\vh 
po Jo V [o,t] 

<CtP{snpE). (6.25) 

[o,t] 

For the third term on the right-hand side of (|6.24p . 

ll:J-[po^T-[C]lliW < ^ll-[/'o*Tlk-ll^"llo 

po V Po 

<C[\\poi"'\\o + \\v\\2]\\ fi"\\ 
Jo 

< CtP{supE) . (6.26) 

[o,t] 



Combining these estimates provides the inequality 



1 / .■! 

supll — [pgw'l |L <CtP{supE) + CsupEi {Mo + tP {sup E)), 
[o,t] Po [0,t] [0,t] [0,t] 

which leads us immediately to: 

supWpov'" + 3pqv"\\o < CtP{snpE)+CsnpEi {Mo + tP{snp E)) . (6.27) 

[0,f] [O.t] [0,t] [0,t] 
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Now, since for any t G [OjTk], the solution v to our parabolic K-problem is in H^{I), we 
infer that pov'" G L'^{I). We can then apply the same integration-by-parts argument as in 
[H to find that 

\\p,v"' + ipWWl - \\Pnv"'\\l + %p',v"\\l + 3 ^ p,p',[\v"\^]' 

= \\p,i"'\\l + %p'y\\l-ij[p,p'^ + p'S]\i"? 

- \\p^v'''\\l + &\\p',i''\\l-ij^p,p'^\v''\^ (6.28) 

Combined with (|6.27p . this yields: 

sup[||pot"''"'l|o + ||Po^^"||o] <CiP(sup£;)-f Csup£;3 (Mq + tP(supi;)) 

[0,t] [0,t] [0,i] [0,t] 

and thus 

sup[||pow"'||o + Wp'ov'h + IIPoS"||o] < Mo + CiP(sup£;) + Csuppi (Afo, 

[0,t] [0,t] [0,t] 

With (O) and (fTSl), it follows that 



sup[||/9ow"'||o + \\i"h] < Mo + CtP{snpE) + C snp (Mo + tP{snpE)), 

[0,t] [0,t] [0,t] [0,t] 

and hence 

sup[||pow"'||o + ll^lb] < A/o + CtP(supi;)-f Csupi;3 (A/q + tP (sup £;)). (6.29) 

[0,t] [0,t] [0,t] [OA] 

6.5. Elliptic and Hardy-type estimates for dfv{t) and dtv{t). We consider the fourth 
time-differentiated version of (|4.1ap : 

which can be rewritten as 

-2 [^^] " 4pldtv']' = -podtv + ci ] + C2 [^^5^] , 

for some constants ci and C2. By employing the fundamental theorem of calculus and 
dividing by pg , we obtain the equation 

rK^tt'J - —[Podtv] =-^dtV + —[ — — J + — [ --,5 J 

Po 
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For any t G [0, T^], Lemma l3^ provides the K-independent estimate 



sup — p^^^V < sup « + sup — ^^^4^4 o+sup — P^T^ 



+ sup||4[Poat'^']'llol|3 / 4lU-+6sup||p|a,V^||„ . (6.30) 

[0,t] p2 7o ^ [0,t] ^ 



We estimate each term on the right-hand side of (|6.30p . The first term on the right-hand 
side is bounded by Mq + Ct P(sup[Q E) thanks to (|6.10p . For the second term on the 
right-hand side of (|6.30p . we have that 



I 1 i Ppdfv V ,n a^-^vu II ^ II I II ^^2~f[ Po V VII 

Po ' ' ' 

< C\\^,{p,d^,i')'Ui'\k + WVp-od^ti'i^] llo + WVp-^djv'iP^] \l 

I .ill ^a2~i\PQ ^'^"l II 

<ciivA^(pof^2)'+ / ^^{p^d^i')'\w\\hv'i+ 1 m 

Jo Ja 2 



"'0 

+ C\\^d^v'\\o\\pov"h+C\\i'h\\ I v"U\pli'2 + I PSd'fi'h 



<\\Vp^{pov'2y+ / Vp^{pov'uyMi'\\!H+ / ^ill! 

Jo Jo = 

r/ iLiir/ii 2 luv _L / ^-/iis 



+ / Vp^Kt\\o\\v\\!H+ / «;||! 

Jo Jo 

+ C\\^d?v'\\o\\pov"\\h\PoVo+ I P^v'IU 





+ C\\i'\\\\\i[+ f i'MW f «"||o|l4^2+ / pldfi'Wi. 
Jo ^ Jo Jo 



where we have again used the fact that || • H^oo < C*!! • ||3. Thanks to the definition of E, 
the previous inequahty shows that for any t G [0,T„], 



sup||^[^2^P^]'|| <C'supi?l {Mo + tP{snpEj). (6.31) 
[o,t] Vpo ^ [o,t] [o,t] 
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For the third term on the right-hand side of ()6.30|) . we have similarly that 

" ' [^T]lW<2|l(poWllo||VA^§flL.+5||v7^^ 



<C\\ipov[y+ f\podttvy\\o\\VP^dti'h + C\\dtv'\\l\\ fipoi'Jh 
Jo Jo 

<c\\ipov[y+ [ ipodui'y\\o\\vp^dti'\\i'''\\i^dti'y\\h-. 

Jo 

+ c\\dti'\\i\\ f\pov"y\\o 



<c\\{poi[y+ l\podui'y\\oH + f dui'wi-m^odtv'ywi.-. 

Jo Jo 

+ C\\dti'\\l\\ [\pov"y\\o, (6.32) 
Jo 

where < a < ^ is given and < a = < 1- 

The only term on the right-hand side of (|6.32[) which is not directly contained in the 
definition of E is \\{y/podtv'y\\22-a- To this end, we notice that 

Wi^dtv'yh^-. < w^h^-^ + wvf^K'tWo 

^VPo 



< \\-^\\^2-i\\dtv'h + \\^i[',\\o, (6.33) 
^VPo 



where we have used the fact that |j • ||lp < Cp || • ||i., for all 1 < p < oo. Thanks to the 
definition of E, the previous inequality and (|6.32[) provides us for any t g [0, T^] with 

sup|| — [4^]'|L < Csup^;" iMo + tP{snpE)), (6.34) 
[o,t] Vo V [o,t] [o,t] 

where we recall again that < a = f^ff < 1- 

The fourth term on the right-hand side of (|6.30p is easily treated: 

^ ' '^*'^']'llo|| f ^\\L^{t) <C[\\pldfv"\\o + \\Vp^d!v'\\o]tsnp\\vh 



Po Jo V [o,t] 

< CtP{snpE) . (6.35) 

Similarly, the fifth term on the right-hand side of (|6.30l) is estimated as follows: 



\Pod'tv'^lit)<C\\p-^dti'\\L^W'\\o 

<C[ \\pld!v"\\o + \\^d!v'\\o]\\fi"\\, 

Jo 

< CtP{supE) . (6.36) 

[o.t] 

Combining the estimates (|6.30p - (|6.36p . we can infer that 

snp\\\[pldfv']'\\ <CtP{supE) + CaupE°' {Mq + tP {sup E)). (6.37) 
[0,t] p2 [o,t] [o,t] [o,t] 
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Now, since for any t e [0, Tk], solutions to our parabolic K-problem have the regularity 
dfV E H^{I), we integrate-by-parts: 



\\^[pld!v']X = \\pld?v''\\l + A\\poh'od?v'\\l + '2 I pUMv^]' 

pI J' 

= \\pldli"\\l + ^p',pld^,v'\\l ^Aj^p'o'poldfif -2j^p'U\d!v'\' 

= \\pld!i'ro-2j^p'U\dM'- (6.38) 

Combined with (|6.37p . and the fact that podfv' = pqv^ + dotpodfv' for the second term 
on the right-hand side of (I6.38|l . we find that 

sup\\pldfi"\\o <CtP{supE) + CsvLpE°' {Mq + tP{sup E)). (6.39) 

[0,t] [O.t] [O.t] [0,t] 



Now, since -^[pldfv']' = p§dfv" + 2p(i^ p'^^d^v' , the estimates (|07l) and dOO]) also imply 

Po 

that 

sup Wpo^p'qO^v'Wo < CtP{supE) + CsupE°' {Mq + tP{supE)). (6.40) 

[0,t] [0,t] [0,t] [0,t] 

Therefore, 

snp[\\pld!v"\\o + Wp'opldfv'Wo + Wpldfi'Wo] <Mo + CtPisnpE) 

[0,t] [O.t] 



CsupE"" {Mo + tP {sup E)), 

[0,t] [0,t] 



so that with ([T3]) and PTB]) . 



anpiWp^dtv'Wo + llpo'^t ^'llo] < Mo + CtP(sup£;) + Csup^;" (Mq + tP(sup£:)) . 

[0,t] [0,t] [0,t] [0,t] 

Together with (|6.13p and the weighted embedding estimate (|2.ip . the above inequality shows 
that 

sup[\\pldfv"\\o + \\dfv\\i] < MQ + CtP{supE) + CsupE'' {Mq + tP{sup E)). (6.41) 

[0,t] ^ [0,t] [0,t] [0,t] 

By studying the second time-differentiated version of (|4.1ap in the same manner, wc find 
that 

sup[||p|9tu"'|lo + WdtvWi] < Mo + CtP{snpE) + CsnpE°' {Mq + tP{siip E)) . (6.42) 

[0,t] ^ [0,t] [0,t] [0,t] 

7. Proof of Theorem If .fl 

7.f. Time of existence and bounds independent of e and existence of solutions 

to (fTTTO]) . Summing the inequalities (jOU)) . ([O^ . (p?^ . (ICTl) . ([O^ . we find that 

sup E{t) < Mo + CtP{aupE) + CsnpE" {Mq + tP{supE)) . 
te[o,T] [o.t] [o.t] [o,t] 

As q; < f , by employing Young's inequality and readjusting the constants, we obtain 

sup E{t) <Mq + CTP{ sup E{t)) . 
te[o,T] te[o,T] 
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Just as in Section 9 of [3j , this provides us with a time of existence Ti independent of k 
and an estimate on (0, Ti) independent of k of the type: 

sup E{t) < 2Mo. (7.1) 

tG[0,Ti] 

In particular, our sequence of solutions (w^) satisfy the K-independent bound (17. ip on the 
K-independent time- interval (0,Ti). 

7.2. The limit as k ^ 0. By the K-independent estimate (|7.1[) . there exists a subsequence 
of {fft} which converges weakly to v in L^(0, T; H^{I)). With ri{t, x) = x + v{s, x)ds, by 
standard compactness arguments, we see that a further subsequence of a-nd ?7j!^ uniformly 
converges to v and r;', respectively, which shows that v is the solution to p.lOp and f (0, x) = 
uo{x). 

7.3. Uniqueness of solutions to the compressible Euler equations (jl.iop . For 

uniqueness, we require the initial data to have one space-derivative better regularity than 
for existence. Given the assumption (jl.lSp on the data (uo,po)j repeating our argument for 
existence, we can produce a solution v on [0, Ti] which satisfies the estimate 



s=0 



< oo . 



and has the flow ri{t, x) — x + v{s, x)ds. For the sake of contradiction, let us assume that 
w is also a solution on [OjTi] with initial data (uo,po), satisfying the same estimate, with 
flow ip{t, x) = X + w{s, x)ds. 
We define 

Sv = V — w , 

in which case we have the following equation for Sv: 

PoSvt + (PoK"' - V^'"'])' - in / X (0, Ti] , (7.2a) 
6v = on / X {< = 0} , (7.2b) 
po = ondl. (7.2c) 

By considering the fifth time-differentiated version of (|7.2ap and taking the L'^{I) inner- 
product with dt5v, we obtain the analogue of (|6.10p (with k — 0) for 5v. The additional 
error terms which arise are easily controlled by the fact that both v and w have one space- 
derivative better regularity than the energy function E. This produces a good bound for 
dfSv G L°°{0,Ti; L^{I)). By repeating the elliptic and Hardy-type estimates for dfSv € 
i:°°(0,Ti;iJi(/)) and dSv € L'^{0,Ti; H^{I)), and using jLlB, we obtain the inequality 

sup {\\dtSv{t)\\l + \\dfSv{tm + \\5v{tm) 
te[o,Ti] 

<CTiP( sup {\\dt5vml + \\d^6v{t)\\l + \\6v{t)\\l)), 
te[o,Ti] 

which shows that Sv = 0. 

7.4. Optimal regularity for initial data. We smoothed our initial data (uqjPo) in or- 
der to construct solutions to our degenerate parabolic K-problem (|4.ip . Having obtained 
solutions which depend only on E(0, v), a standard density argument shows that the initial 
data needs only to satisfy Mq < oo. 
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8. The case 7 7^ 2 

In this section, we describe the modifications to the energy function and the methodology 
for the case that 7^2. We denote by qq the integer satisfying the inequaUty 

1< 1 + - - ao < 2 . 

7- 1 

Letting 

d{x) = dist(a;, dl) , 
We consider the foUowing higher-order energy function: 
4 2 

E,{t, v)^J2 •)ll2-f + E \\dd?'^it, Ollts + \\Vdd,dlv{t, + WVdd^dMt, Olio 

s=0 s=0 
ao 



|2 



' 

and define the polynomial function = P{Ej{0,v)). Notice the last sum in E-y appears 
whenever 7 < 2, and the number of time-differentiated problems increases as 7 — > 1. 

Using the same procedure as we have detailed for the case that 7 = 2, we have the 
following 

Theorem 8.1 (Existence and uniqueness for any 7 > 1). Given initial data {uo,Pq) such 
that Mq < 00 and the physical vacuum condition Jj.^p holds for pQ, there exists a solution 
to lll.9\) (and hence on [0,T-y] for > taken sufficiently small, such that 

sup E{t) < 2M^ . 
te[o,T] 

Moreover if the initial data satisfies 

3 3 ao J 

E 11^*^(0' ■)\\H-^-^ii} + E WddXO, Wm-^ii) + E 11^ ^-^^'^9r''«-V(0, Oll^ < 00 , 

s=0 s=0 a=0 

then the solution is unique. 
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